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Optimization  of  composite  structures  subject  to  local  stress 

constraints 


Robert  Lipton  ^and  Michael  Stuebner 
Mathematics  Department,  Louisiana  State  University, 
Baton  Rouge,  Louisiana  70803 


Abstract.  An  extension  of  current  methodologies  is  introduced  for  optimization  of 
graded  microstructure  subject  to  local  stress  criteria.  The  method  is  based  on  new  multi¬ 
scale  stress  criteria  given  by  macrostress  modulation  functions.  The  modulation  functions 
quantify  the  intensity  of  local  stress  fluctuations  at  the  scale  of  the  microstructure  due  to 
the  imposed  macroscopic  stress.  The  methodology  is  illustrated  for  long  cylindrical  shafts 
reinforced  with  stiff  cylindrical  elastic  fibers  with  generators  parallel  to  the  shaft.  Exam¬ 
ples  are  presented  for  shaft  cross  sections  that  possess  reentrant  corners  typically  seen  in 
lap  joints  and  junctions  of  struts.  It  is  shown  that  the  computational  methodology  delivers 
graded  fiber  microgeometries  that  provide  overall  structural  rigidity  while  at  the  same  time 
tempering  the  influence  of  stress  concentrations  near  reentrant  corners. 

Keywords:  Stress  constraints,  optimal  structural  design 


1  Introduction 

The  high  specific  stiffness  of  composite  materials  makes  them  attractive  structural  materials 
for  use  in  aerospace  applications.  These  materials  are  increasingly  being  used  in  geometries 
that  involve  abrupt  dimensional  changes  within  structural  components,  such  as  skins  con¬ 
nected  to  ribs,  panel  reinforcements  and  junctions  of  struts.  These  geometries  contain  stress 
concentrations  and  thus  present  potential  sources  for  structural  failure.  In  this  paper  we 
present  a  computational  design  method  for  grading  the  microstructure  in  order  to  control 
local  stress  in  the  vicinity  of  stress  concentrations.  The  method  is  based  upon  new  rigorous 
multiscale  stress  criteria  that  connect  the  macroscopic  or  homogenized  stress  to  local  stress 
fluctuations  at  the  scale  of  the  microstructure  [25,  26,  27]. 

It  is  now  well  known  that  effective  constitutive  properties  relating  average  stress  to  av¬ 
erage  strain  can  be  employed  in  the  numerical  design  of  composite  structures  for  optimal 
structural  compliance  and  natural  frequency.  This  type  of  design  problem  has  received 
significant  attention  from  both  the  applied  mathematics  and  structural  optimization  com¬ 
munities  in  the  1980s  and  1990s  [1,  4,  7,  10,  11,  12,  14,  17,  18,  31,  32,  35].  This  list  of 
references  is  by  no  means  complete  and  the  reader  is  referred  to  the  reviews  and  manu¬ 
scripts  [2,  5,  8,  9,  19,  30,  33,  37].  The  associated  body  of  methods  are  referred  to  as  the 
homogenization  approach  to  structural  optimization.  These  treatments  have  focused  on  the 
design  problem  in  the  absence  of  stress  or  strain  constraints. 

'^Corresponding  Author:  E-mail  address:  lipton@math.lsu.edu,  Phone:-|- 1-6176-538782,  Fax:-|-l-2255- 
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Recent  efforts  have  initiated  the  development  of  numerical  methods  for  structural  op¬ 
timization  in  the  presence  of  stress  constraints.  The  investigation  given  in  [13]  provides  a 
numerical  method  for  the  stress  constrained  minimum  volume  design  problem.  The  method 
is  carried  out  using  an  empirical  model  that  is  an  extension  of  the  power  penalized  stiffness 
model  also  known  as  the  Solid  Isotropic  Microstructure  with  Penalization  (SIMP)  model 
[5].  The  choice  of  local  stress  constraints  proposed  in  [13]  is  motivated  by  the  explicit 
form  of  the  corrector  tensors  associated  with  rank  two  orthogonal  laminar  microstructures. 
Other  work  has  extended  the  numerical  methodology  providing  the  opportunity  for  design 
of  composites  and  topology  in  the  presence  of  mean  square  stress  constraints.  The  prob¬ 
lem  of  mean  square  stress  constrained  structural  optimization  for  fiber  reinforced  shafts  is 
taken  up  in  [20].  In  that  work  a  numerical  algorithm  is  developed  based  on  a  suitable  ho¬ 
mogenized  quantity  (the  covariance  tensor)  that  accurately  encodes  the  mean  square  stress 
constraints.  The  work  of  [3]  introduces  a  partial  relaxation  for  topology  optimization  for 
minimum  mean  square  stress  using  finite  rank  laminates.  The  rigorous  theoretical  con¬ 
text  behind  the  homogenization  approach  to  mean  square  stress  (or  gradient)  constrained 
structural  optimization  has  been  worked  out  in  [23],  [24].  Other  recent  theoretical  and  nu¬ 
merical  developments  have  addressed  problems  of  minimizing  objective  functions  depending 
explicitly  on  the  mean  square  norm  of  the  gradients  of  the  state  variable  [15,  22,  34,  36,  38]. 

In  this  treatment  we  consider  design  problems  that  require  point  wise  constraints  on 
the  local  stress.  Suitable  multiscale  quantities  beyond  effective  constitutive  laws  are  used 
to  accurately  enforce  local  stress  criteria.  These  quantities  dubbed  macrostress  modulation 
functions  are  introduced  in  [25,  26,  27].  The  design  formulation  given  here  is  expressed  in 
terms  of  a  homogenized  design  problem  that  satisfies  two  requirements:  The  first  is  that  the 
homogenized  design  problem  is  computationally  tractable.  The  second  is  that  the  solution 
of  the  homogenized  design  problem  provides  the  means  to  identify  graded  microstructures 
that  deliver  the  required  structural  response  while  at  the  same  time  provide  local  stress 
control.  The  approach  put  forth  in  this  paper  is  in  essence  an  inverse  method  that  uses  new 
tools  from  homogenization  theory  to  identify  optimal  graded  microgeometries.  Because  of 
this  we  will  refer  to  it  as  an  inverse  homogenization  design  method.  The  rigorous  theoretical 
basis  for  this  approach  is  given  in  [21]  and  [29]. 

To  illustrate  the  ideas  this  article  treats  the  problem  of  reinforcement  of  a  long  shaft  with 
constant  cross  section  subjected  to  torsion  loading.  The  microstructure  within  the  shaft 
consists  of  long  reinforcement  fibers  of  constant  cross  section  with  isotropic  shear  modulus 
Gf  embedded  in  a  more  compliant  material  with  shear  modulus  Gm-  The  shaft  together 
with  the  fibers  are  right  cylinders  with  generators  along  the  X3  axis.  The  cross  section  of 
the  reinforced  shaft  is  specified  by  the  region  0,  in  the  xi  —  X2  plane.  In  the  neighborhood  of 
any  point  x  =  (xi,  X2)  the  local  microgeometry  is  given  by  a  periodic  geometry  with  period 
cell  filled  with  the  Hashin  Shtrikman  coated  cylinder  assemblage  [16],  see  Figure  1.  The 
characteristic  length  scale  of  the  period  relative  to  the  size  of  the  design  domain  is  taken 
to  be  small  and  is  denoted  by  e.  The  coated  cylinder  assemblage  corresponds  to  stiff  fibers 
of  circular  cross  section  separated  by  compliant  matrix  material.  The  local  area  fraction  of 
the  fiber  phase  is  given  hy  0f  and  the  matrix  area  fraction  is  given  by  9m  with  9f  +  9m  =  1  • 
The  graded  material  properties  of  the  shaft  are  obtained  by  changing  the  local  area  fraction 
of  fibers  across  the  cross  section  of  the  shaft  see  Figure  1. 

In  what  follows  a  constraint  is  placed  on  the  total  cross-sectional  area  occupied  by  the 
fibers.  The  goal  of  the  design  problem  is  to  identify  a  graded  distribution  of  fibers  across 
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the  cross  section  such  that  the  following  requirements  are  met: 

I.  The  reinforced  shaft  has  a  torsional  rigidity  that  is  acceptable. 

II.  The  magnitude  of  the  local  point  wise  stress  inside  the  composite  is  controlled  over  a 
designated  subset  of  the  cross  section. 

In  Section  2  we  describe  a  suitably  formulated  homogenized  design  problem  that  allows  for 
the  easy  identification  of  fiber  reinforced  designs  satisfying  requirements  (I)  and  (II).  The 
homogenized  formulation  makes  use  of  a  multiscale  stress  criterion  given  in  terms  of  the 
macrostress  modulation  function.  The  physical  motivation  for  the  multiscale  stress  criterion 
is  presented  in  Section  2.  The  local  effective  properties  and  macrostress  modulation  func¬ 
tion  used  to  formulate  the  stress  criterion  are  computed  explicitly  for  the  Hashin  Shtrikman 
coated  cylinder  assemblage  in  [26].  These  are  given  in  terms  of  the  elastic  moduli  of  the 
matrix  and  fibers  and  the  local  area  fraction  of  fibers.  We  point  out  that  the  coated  cylin¬ 
der  assemblage  gives  the  minimum  amount  of  stress  amplification  among  all  periodic  fiber 
microstructures  with  isotropic  effective  shear  compliance  [28].  A  gradient  based  numerical 
implementation  of  the  inverse  homogenization  design  method  is  described  in  Section  3.  In 
Section  4  computational  results  are  given  for  the  “X”  shaped  shaft  cross  section.  This  geom¬ 
etry  typifies  the  junctions  between  composite  substructures  and  possesses  reentrant  corners 
seen  in  lap  joints  and  junctions  of  struts.  We  conclude  by  noting  that  the  theoretical  basis 
for  the  approach  given  here  has  been  established  for  three  dimensional  structural  design 
using  multiphase  locally  periodic  composites  in  the  presence  of  point  wise  stress  constraints 
see  ([21],  Theorems  5.1  and  5.2).  For  locally  layered  microstructures  the  corresponding 
theory  is  presented  in  [29] . 

2  Homogenized  design  formulation  and  identification  of  op¬ 
timal  graded  fiber  microgeometries. 

The  inverse  homogenization  design  method  is  a  top  down  design  approach.  First  a  well 
posed  homogenized  design  problem  is  developed.  This  design  problem  is  given  in  terms  of 
design  variables  that  reflect  the  local  microgeometry  inside  the  composite.  For  the  problem 
treated  here  the  design  variable  for  the  homogenized  design  problem  is  given  by  the  density 
function  0/(x).  The  homogenized  design  problem  is  then  solved  to  obtain  an  optimal 
density  function.  With  the  optimal  density  in  hand  we  use  it  to  recover  an  explicit  graded 
fiber  design  that  has  structural  properties  close  to  that  of  the  optimal  homogenized  design 
and  satisfies  prescribed  point  wise  stress  constraints.  The  homogenized  design  problem  is 
described  in  the  first  subsection.  The  explicit  link  between  homogenized  designs  and  graded 
fiber  reinforced  designs  that  satisfy  point  wise  stress  constraints  is  provided  in  the  second 
subsection. 

2.1  Homogenized  design  problem 

The  design  variable  for  the  homogenized  design  problem  is  given  by  the  density  function 
0j(x).  This  function  is  interpreted  as  providing  the  local  area  fraction  of  the  fiber  phase  in 
a  homogenized  design.  The  resource  constraint  on  the  fiber  phase  is  given  by 

/  9 f{'x)  dxidx2  <  0  X  {Area  of  H),  (2-1) 

Jn 
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where  0  <  0  <  1.  At  each  point  the  local  area  fraction  satisfies  the  box  constraint  given  by 

0<9f<l.  (2.2) 

In  this  treatment  the  local  fiber  area  fraction  6 f  changes  continuously  with  position  accord¬ 
ing  to  the  condition 


|6»/(x)  -  6i/(x-Fh)|  <  A:|h|. 


(2.3) 


Here  the  constant  K  is  prescribed  by  the  designer.  The  universe  of  admissible  designs  given 
by  all  local  area  fractions  9f  satisfying  the  resource  constraint,  box  constraints,  and  (2.3) 
is  denoted  by  Dq. 

The  compliance  in  shear  for  the  matrix  and  fiber  are  given  by  Sm  =  (2Gm)“^  and 
Sf  =  {2Gf)~^  respectively.  Here  the  matrix  is  more  compliant  and  Sm  >  Sf.  The  effective 
shear  compliance  S^{9f)  for  the  Hashin  Shtrikman  coated  sphere  assemblage  made  from 
stiff  fibers  with  area  fraction  9f  is  given  by  [16] 


5^(%(x)) 


/Sm  +  Sf+9f{^){Sf-Sm)\ 
\Sm  +  Sf+9f{^){Sf  +  Sm))  ■ 


(2.4) 


The  macroscopic  stress  potential  vanishes  on  the  boundary  of  the  shaft  cross  section 
and  satisfies 

-div  {S^{9f)V<j)^)  =  1  (2.5) 

inside  the  cross  section.  The  torsional  rigidity  for  the  homogenized  shaft  cross  section  made 
from  a  homogenized  material  with  compliance  S^{9f)  is  given  by 


TZ{9f)  =  2  I  dxidx2- 


IQ. 


(2.6) 


The  macroscopic  stress  in  the  homogenized  shaft  is  given  by  =  KVcj)^  where  R  is  the 
rotation  matrix  associated  with  a  counter  clock  wise  rotation  of  7r/2  radians. 

The  multiscale  stress  criterion  is  given  in  terms  of  the  macrostress  modulation  function 
introduced  in  [25].  The  macrostress  modulation  function  captures  the  interaction  between 
the  macroscopic  stress  cj^(x)  and  the  microstructure.  Consider  the  unit  square  Q  filled 
with  Hashin  Shtrikman  coated  cylinder  assemblage  with  area  fraction  of  fibers  9f.  The 
coordinates  of  points  inside  the  unit  square  are  denoted  by  y.  The  associated  local  shear 
compliance  inside  the  unit  square  is  denoted  by  S{9f,y).  Here  S{9f,y)  =  Sf  in  the  fiber 
phase  and  S{9f,y)  =  Sm  in  the  matrix  phase.  In  what  follows  all  derivatives  with  respect 
to  the  microscopic  variable  y  are  denoted  with  subscripts.  The  microscopic  response  to  the 
imposed  macroscopic  stress  is  given  by  cj(x,  y)  =  i?[Vy(rc(x, y))  -|- V(^^(x)],  where  for  each 
X  in  the  shaft  cross  section  H  the  Q  periodic  fluctuating  stress  potential  rc(x,  y)  solves  the 
microscopic  equilibrium  equation 


-divy  (5(6'/(x),y)(Vy(u;(x,y)) -h  V(()^(x)))  =  0,  y  in  Q.  (2.7) 


Here  the  x  coordinate  appears  as  a  parameter.  The  relevant  interaction  is  described  by  the 
macrostress  modulation  function  f{9f,a^)  given  by 

/((9/(x),fT^(x))  =  sup  {lfT(x,y)l2}  .  (2.8) 

y  in  Q 


4 


Physically  the  macrostress  modulation  provides  an  upper  envelope  on  the  oscillating  point 
wise  local  stress  in  the  composite  [25] ,  [26] . 

The  macrostress  modulation  is  calculated  explicitly  for  the  Hashin  Shtrikman  coated 
cylinder  assemblage  in  [26].  We  define  the  amplification  factor 

^  {sm  +  Sf  +  ef{Sm-Sf)) 

and  set 


/(6*/,v)  =  A{9f)\\-\‘^  if  6f  >  0  and  =  Jvj^,  if  =  0  (2-10) 

for  every  vector  v.  Here  A{9f)  >  1  and  H(0)  =  {2Sm/{Sm  +  Sf))"^  >  1.  In  the  context  of 
torsional  rigidity  the  macrostress  modulation  is  written  in  terms  of  the  homogenized  stress 
potential  and  is  given  by 


/(6»/(x),  V(/)^(x)).  (2.11) 

It  is  pointed  out  that  the  Hashin  Shtrikman  coated  cylinder  assemblage  gives  the  smallest 
amplification  factor  A{9f)  among  all  locally  periodic  fiber  microstructures  with  isotropic 
effective  shear  compliance  [28]. 

We  choose  a  subset  S  of  the  shaft  cross  section  that  lies  a  finite  distance  away  from  the 
boundary.  On  this  set  the  prescribed  multiscale  stress  criterion  is  given  by 

/(0/(x),V<()^(x))  <r2.  (2.12) 

In  this  treatment  domains  with  reentrant  corners  are  considered  so  there  will  be  a  stress 
singularity  at  the  corner.  Therefore  the  choice  of  T  >  0  depends  on  the  distance  between  S 
and  the  reentrant  corner.  It  is  clear  that  the  multiscale  stress  criterion  may  not  be  satisfied 
by  any  homogenized  design  if  T  is  chosen  too  small. 

The  homogenized  design  problem  is  given  by 

HP  in!  ms,)}-,  subject  to:  |  tor  x  in  5  }  }' 

The  homogenized  design  problem  HP  is  well  posed  and  there  is  an  optimal  design  9f 
provided  T  is  chosen  large  enough  so  that  Dq  contains  at  least  one  design,  see  [21]. 

2.2  Identification  of  graded  fiber  design  from  the  homogenized  design 

In  this  subsection  it  is  shown  how  to  use  the  optimal  design  9 j  of  HP  to  identify  a  graded 
fiber  design  satisfying  the  requirements  (I)  and  (H).  The  building  block  for  the  microstruc¬ 
ture  is  the  unit  cell  filled  with  a  Hashin  Shtrikman  coated  cylinder  assemblage  [16].  The 
coated  cylinder  assemblage  is  constructed  as  follows.  A  space  filling  configuration  of  disks 
of  different  sizes  ranging  down  to  the  infinitesimal  is  placed  inside  the  unit  period  cell  given 
by  the  unit  square  Q.  Each  disk  is  then  partitioned  into  an  annulus  called  the  coating  and 
a  concentric  disk  which  is  a  fiber  cross  section.  The  union  of  these  “coated  disks”  make 
up  the  coated  cylinder  assemblage.  The  area  fraction  of  the  fiber  phase  is  the  same  for  all 
coated  disks  in  the  assemblage  and  is  given  hy  9f.  The  union  of  the  coatings  comprises  the 
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matrix  phase.  The  area  fraction  of  the  fiber  phase  for  the  coated  cylinder  assemblage  inside 
the  unit  square  Q  is  easily  seen  to  be  Of.  The  unit  square  filled  with  the  coated  cylinder 
assemblage  is  illustrated  in  Figure  1.  A  microstructure  is  obtained  by  rescaling  the  unit  cell 
by  the  factor  e  so  that  it  becomes  the  period  cell  for  an  e  periodic  composite. 

In  order  to  describe  the  graded  fiber  composite  one  partitions  the  shaft  cross  section 
into  the  N  subdomains  k  =  1, ...  ,N  and  to^.  The  maximum  diameter  of 

the  subdomains  in  the  partition  is  denoted  by  .  The  partition  is  denoted  by  .  A 
graded  fiber  composite  is  constructed  by  placing  an  e  periodic  Hashin  Shtrikman  coated 
cylinder  geometry  inside  each  subdomain  uj^.  The  area  fraction  of  fibers  in  each  subdomain 
is  given  by  the  constant  0^  and  these  constants  can  change  between  subdomains.  For  future 
reference  this  type  of  locally  periodic  microstructure  will  be  called  a  (tat,  e)-graded  Hashin 
Shtrikman  fiber  microstructure. 

The  local  piecewise  constant  shear  compliance  for  the  (tat,  e)-graded  Hashin  Shtrikman 
fiber  microstructure  is  denoted  by  .  The  stress  potential  for  this  microstructure  is 
denoted  by  and  vanishes  on  the  boundary  of  the  cross  section.  The  stress  potential 
satisfies  the  equilibrium  equation 

-div  =  1.  (2.14) 

The  torsional  rigidity  of  the  cross  section  is  given  by 

7^^A  =  2/'  dxidx2.  (2.15) 

Jn 

The  nonzero  components  of  the  in  plane  stress  are  denoted  by  the  vector  = 

,  o'23^)  and  are  related  to  the  gradient  of  the  stress  potential  according  to 

=  RVcj)^’^.  (2.16) 

Here  R  is  the  matrix  corresponding  to  a  counter  clockwise  rotation  of  tt/2  and  = 

Given  any  partition  the  partition  is  said  to  be  a  refinement  of  it,  if  <  tm 
and  if  each  subdomain  belonging  to  is  a  subset  of  a  subdomain  in  .  Next  for 
any  given  partition  we  introduce  a  nested  sequence  of  refinements  such 

that  and  is  a  refinement  of  with  lim^^oo  =  0.  Here  we 

assume  e  <  Tat^  and  we  will  consider  the  sequence  of  (rAr<;, e)-graded  Hashin  Shtrikman 
fiber  microstructures  with  associated  stress  potentials  (x). 

The  identification  of  a  graded  fiber  composite  satisfying  the  point  wise  stress  constraints 
is  given  in  the  following  Proposition. 


Proposition  2.1.  Identification  of  graded  microstructure. 

Consider  a  homogenized  design  speeified  by  Of  in  Dq  sueh  that  the  multiseale  stress 
eriterion  (2.12)  holds  on  the  subset  S  of  the  shaft  eross  seetion.  Consider  also  any  partition 
of  the  design  domain  H.  Then  for  any  given  t  >  T  and  small  number  (5  >  0,  there  is  a 
refinement  of  the  partition  and  an  assoeiated  graded  Hashin  Shtrikman  fiber 

mierostrueture  for  whieh  the  part  of  S  over  whieh  the  stress  eonstraint 

|V((."’^'(x)|  <t  (2.17) 
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is  violated  has  measure  (area)  less  than  5  and 

-7^(0/)|  <  5,  (2.18) 

and 

Nt 

<e  X  {Area  of  n)  +  6.  (2.19) 

k=l 

Here  \u!^\  denotes  the  area  of  .  Moreover  the  area  fraetions  9^  are  determined  from  Of 
through  the  averages  given  by 

9f^  =  X  f  9f{x)dxidx2-  (2.20) 

This  Proposition  is  established  in  [21].  The  homogenized  design  formulation  together 
with  Proposition  2.1  provide  an  inverse  homogenization  method  for  identifying  microstruc¬ 
tures  that  satisfy  point  wise  stress  constraints  while  delivering  a  torsional  rigidity  close  to 
that  given  by  the  optimal  design  9f  for  the  homogenized  design  problem. 

3  Gradient  algorithm  for  the  homogenized  design  problem 

In  the  computational  examples  we  enforce  the  stress  constraint  by  adding  a  penalty  term 
to  the  torsional  rigidity  and  minimize 

LiOf)  =  -n{9f)  +  I  [  if  {Of,  dxidx2,  (3.1) 

Jn 

over  all  Of  in  Dq  where  I  >  0  and  satisfies 

-div  {S^{0f)V(j)^)  =  1  (3.2) 

and  vanishes  at  the  boundary.  The  computational  examples  provided  here  will  be  carried 
out  for  a  domain  with  reentrant  corners  of  interior  angle  Sir (2.  In  view  of  the  strength  of 
the  associated  singularity  at  the  reentrant  corners  the  power  “p”  appearing  in  the  penalty 
term  is  chosen  to  be  less  than  3.  This  minimization  problem  is  well  posed  and  one  can  find 
a  minimizing  density  0/(x)  see  [21]. 

Given  the  minimizing  density  and  associated  stress  potential  (j)^  one  considers  sets  of 
the  form 

=  {x  G  P!  :  f{9f{x),V^^{x))  <  T2}  (3.3) 

For  fixed  choices  of  d  >  0  and  t  >  T  Proposition  2.1  can  be  used  to  recover  a  (r^,  e)-graded 
Hashin  Shtrikman  fiber  microstructure  for  which  the  the  part  of  At  over  which  the  stress 
constraint 

|V</.^’^(x)l<t  (3.4) 

is  violated  has  measure  less  than  6, 

_7^(0^)|  <  (3.5) 
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and 


N 

<e  X  {Area  of  n)  +  S.  (3.6) 

k=l 

Here  the  area  fractions  9^  are  given  by 

=  -r^  X  f  9f{x)dxidx2-  (3.7) 

1^  I  J  UJ^ 

To  proceed  numerically  we  compute  the  derivative  of  the  objective  function  given  by 
(3.1)  with  respect  to  the  design  variable  Of.  The  derivative  will  be  used  to  develop  a  gradient 
minimization  algorithm.  However  we  note  that  the  macrostress  modulation  function  (2.10) 
is  discontinuous  at  6*/  =  0.  In  order  to  proceed  we  restrict  Of  to  satisfy  the  constraint 
0  <  0™*”'  <  Of  <  1  at  each  point  in  the  design  domain.  For  our  computations  we  choose 
Qjiin  _  Q  Q]^_  Pqj.  choice  the  modulation  function  is  differentiable  with  respect  to  the 
density  and  we  proceed  to  develop  a  gradient  method  subject  to  these  new  constraints. 

The  gradient  of  the  objective  is  computed  with  the  help  of  an  adjoint  held  A.  Here  A  is 
the  solution  of 


-div  (5®(0/) VA)  =  1  +  /  div  {2p{A{0 f)V ■  V(j)^)P-^A{0f)V(p^)  .  (3.8) 

where  1  <  p  <  3  and  A  =  0  on  the  boundary.  For  77  <<  1  the  change  in  the  stress  potential 
due  to  small  local  perturbations  rjOf  in  the  area  fraction,  is  written  as  (p  and 

-div  (^5^(0/)V(^)  =  div  (^{d0^S^{Of)Of)V(p^^  (3.9) 

where  0  =  0  on  the  boundary.  Applying  standard  arguments  using  (3.2),  (3.8)  and  (3.9) 
one  calculates  to  hrst  order  that 


AL 

dojL 


rj  /  de.LOf  dxidx2, 


where 


dejS^{0f)VX  ■  V(p^  +  2lp{A{ef)V(p^  ■  V(p^)P-^defA{0f)V<p^  ■  V(p^{3.10) 


The  continuity  constraints  on  0f{x)  expressed  by  (2.3)  are  enforced  by  the  way  in  which 
the  design  variable  is  initialized  and  updated.  The  local  average  of  a  scalar  function  /  over 
the  disk  of  radius  R  centered  at  p  is  denoted  by  <  /  (p).  For  a  given  held  Of  satisfying 

the  resource  constraint  (2.1)  and  box  constraint  0™*"'  <  Of  <  1  the  initial  choice  of  design 
variable  9^  is  given  by 


0j  =  {0ff{x). 


(3.11) 


At  the  nth  step  we  suppose  that  is  given  and  we  solve  for  cp^  and  A  using  the  system  of 
equations  (3.2)  and  (3.8).  The  design  variable  0{f  is  updated  according  to 


e]+^  =  {0]-7jde^L)^{x)  (3.12) 

were  de^L  is  given  by  (3.10).  If  the  right  hand  side  of  (3.12)  lies  out  side  the  box  constraint 
gmin  <  <  1  -vve  Update  according  to 


on+i  ^  {0] -rjde^L)^{x)  and 

9]+^  =  1,  if  {0]-pde,L)^{x)>l. 
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(3.13) 


Because  the  updated  functions  are  given  by  averages  of  bounded  functions  it  is  easily  seen 
that  they  satisfy  (2.3)  for  a  non-negative  constant  K  independent  of  x.  The  use  of  local 
averaging  in  the  update  scheme  is  similar  to  the  use  of  filters  in  topology  optimization  see 
[5]  and  [6]. 

For  points  near  the  boundary  a  difficulty  arises  when  defining  the  averages.  This  is  dealt 
with  by  extending  Of  to  the  slightly  larger  domain  17^  =  {x  in  R^;  (iist(x,  17)  <  R}.  The 
particular  form  of  extension  is  up  to  the  designer.  Possibilities  include  setting  Of  =  1  m 
17 \  17  or  reflection  of  Of  across  the  boundary  of  17  into  17 /j.  In  the  discretized  problem 
used  for  the  simulations  we  allow  Of  to  take  constant  values  inside  each  element  and  define 
{Of)^  to  be  the  average  of  Of  taken  over  all  neighboring  elements. 

4  Numerical  results  for  the  X-shaped  cross  section 

The  computational  examples  are  carried  out  for  an  “X”  shaped  domain.  All  interior  angles 
for  the  reentrant  corners  are  fixed  at  3tt/2  radians.  The  shear  stiffness  of  the  matrix  is 
assigned  the  value  Gm  =  1  GPa  and  the  shear  stiffness  of  the  fiber  phase  is  assigned  the 
value  Gf  =  2  GPa.  For  these  choices  Sm  =  l/{2Gm)  =  0.5  and  Sf  =  1/(2G/)  =  0.25. 
All  of  the  design  optimizations  are  carried  out  with  the  matrix  material  occupying  70%  of 
the  shaft  cross  section.  The  local  fiber  area  fraction  Of  is  constrained  to  lie  in  the  interval 

0.01  <0f  <1. 

In  the  first  example  we  optimize  for  torsional  rigidity  only.  The  resulting  optimal  design 
is  referred  to  as  design  1.  The  grey  scale  plot  of  the  local  area  fraction  of  matrix  material 
Om  =  i-—Sf  is  given  in  Figure  2.  Here  the  lightest  regions  correspond  to  points  where  Om  =  0. 
The  darkest  regions  correspond  to  points  where  Om  =  0.99.  As  expected  this  design  ignores 
the  stress  concentration  at  the  reentrant  corners  and  a  zone  of  stiff  material  surrounds  a 
more  compliant  inner  core  and  lies  adjacent  to  the  reentrant  corners.  In  the  next  example 
the  torsional  rigidity  is  optimized  in  the  presence  of  an  integral  penalization  f  (f)^,  i.e., 
p  =  1  for  the  Lagrangian  in  (3.1).  The  resulting  design  is  referred  to  as  design  2.  The  plot 
of  Om  for  this  design  is  given  in  Figure  3.  For  this  case  it  is  seen  that  there  is  a  topology 
change  in  the  design  and  the  more  compliant  material  now  surrounds  the  stiffer  material 
shielding  it  from  the  stress  concentration  at  the  reentrant  corners.  In  the  final  example  the 
torsional  rigidity  is  optimized  in  the  presence  of  a  stronger  integral  penalization  given  by 
f  (/)^,  i.e.,  p  =  2  for  the  Lagrangian  in  (3.1).  The  resulting  design  is  referred  to  as  design 
3.  The  plot  of  Om  for  this  design  is  given  in  Figure  4.  As  in  design  2  it  is  seen  that  the  more 
compliant  material  shields  the  stiff  material  from  the  stress  concentration  at  the  reentrant 
corners.  The  associated  torsional  rigidities  for  each  of  these  cases  are  listed  in  Table  1.  It 
is  seen  from  the  table  that  the  torsional  rigidity  drops  by  44%  for  the  stress  constrained 
designs. 

The  contour  plot  of  the  macrostress  modulation  function  /  for  design  1  is  given  in 
Figure  5.  Figures  6  and  7  give  the  contour  plots  for  /  in  designs  2  and  3  respectively.  When 
comparing  designs  1,  2  and  3  it  is  clear  from  Figures  5,  6  and  7  that  designs  2  and  3  provide 
a  significant  reduction  in  the  size  of  the  over  stressed  zone  /  >  1. 

It  is  pointed  out  that  Proposition  2.1  provides  the  method  for  constructing  a  {tn,s)- 
graded  Hashin  Shtrikman  fiber  microstructure  from  the  data  given  in  design  2.  To  fix 
ideas  we  choose  a  tolerance  5  =  1/1000  and  t  =  1.001.  Then  Proposition  2.1  shows  how 
to  construct  a  (tat,  e)-graded  Hashin  Shtrikman  fiber  microstructure  with  torsional  rigidity 
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satisfying 


-  0.2041  <  1/1000  (4.1) 

and  for  which  the  magnitude  of  the  in  plane  stress  lies  below  1.001  for  all  points  in  the 
region  /  <  1  of  Figure  6,  with  the  possible  exception  of  a  subset  of  points  of  area  less  than 
1/1000. 
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Design  ^ 

Stress  Constraint 

Matrix- Volume  Fraction 

Torsional  Rigidity 

1 

None 

70% 

0.368 

2 

lUY 

70% 

0.204 

3 

jvr 

70% 

0.204 

Table  1: 
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Figure  1:  Cross  section  of  shaft  filled  with  a  graded  locally  periodic  Hashin  Shtrikman 
coated  cylinder  assemblage.  The  change  in  the  local  fiber  area  fraction  across  the  structure 
is  illustrated. 
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The  Final  Density  Distribution 


Figure  2:  Design  1.  Grey  scale  plot  of  the  area  fraction  of  matrix  material  inside  the  X- 
shaped  shaft  cross  section  optimized  for  torsional  rigidity  only.  Lightest  regions  are  made 
exclusively  of  stiff  fiber-phase  material.  Darkest  regions  contain  only  1%  by  area  fraction 
of  fiber  phase  material. 


The  Final  Density  Distribution 
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Figure  3:  Design  2.  Area  fraction  distribution  of  matrix  material  inside  the  X-shaped  cross 
section  optimized  for  torsional  rigidity  with  p  =  1  integral  penalty  on  /. 
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The  Final  Density  Distribution 


Figure  4:  Design  3.  Area  fraction  distribution  of  matrix  material  inside  the  X-shaped  cross 
section  optimized  for  torsional  rigidity  with  p  =  2  integral  penalty  on  /. 


Figure  5:  Contour  plot  of  /  for  design  1. 
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Figure  6:  Contour  plot  of  /  for  design  2. 


Figure  7:  Contour  plot  of  /  for  design  3. 
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